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1 Introduction. 



Interpolating gauges have been proven to be a powerful tool to achieve a better understanding 
of the behavior of the Green's functions in different gauges. Let us quote, for instance, the 
interpolating gauge introduced in [Ij, allowing one to connect the Landau and the light cone 
gauges, thereby providing a useful framework to investigate Yang-Mills theories in noncovari- 
ant gauges. Another example is provided by the interpolating gauge proposed in [21 IS]; which 
connects the Landau and the Coulomb gauges. More recently, a class of covariant gauges in- 
terpolating between the Landau and the maximal Abelian gauge has been introduced in 4 in 
order to investigate the properties of the vacuum energy of Yang- Mills theories due to dimension 

two gauge condensates. More precisely, it has been possible to show that the vacuum energy 

A a A a 

obtained in the maximal Abelian gauge, due to the condensate ( M 2 M +ac a c a ), is the same as the 

A A A A 

vacuum energy obtained in the Landau gauge, due to the condensate* ( M 2 M ) • We recall that 

A a A a 

a nonvanishing condensate (— £ y- £i + ac a c a ) in the maximal Abelian gauge results in a dynam- 
ical mass generation for the off-diagonal gluons. This can be seen as evidence for the Abelian 
dominance hypothesis, implying that Yang-Mills theories in the low energy region should be 
described by an effective abelian theory, according to the dual superconductivity mechanism for 
color confinement. 

In this work we introduce a class of covariant gauges allowing one to interpolate between the 
Landau gauge, the maximal Abelian gauge, the linear covariant gauges and the Curci-Ferrari 
gauge. It is worth remarking that in all these gauges a nonvanishing dimension two condensate 
has emerged (see |U0I] an d references therein), so that this generalized interpolating gauge might 
be useful for further investigation of these condensates. 

The work is organized as follows. In Sect. 2 we introduce the interpolating gauge fixing. 
Since the maximal Abelian gauge requires a splitting of the SU(N) color index A = (i,a), 
corresponding to the N — 1 diagonal generators of the Cartan subalgebra of SU(N), labeled 
by the index i, and to the N(N — 1) off-diagonal generators, labeled by a, it turns out that 
the gauge-fixing term contains six independent gauge parameters. In Sect. 3 we discuss the 
various limits for those gauge parameters, allowing one to recover the linear covariant gauges, 
the Landau gauge, the maximal Abelian gauge and the Curci-Ferrari gauge. In Sect. 4 we 
prove the all orders multiplicative renormalizability by means of the algebraic renormalization 
procedure. In Sect. 5 we present the explicit one-loop calculation of the anomalous dimensions 
of all fields and parameters. These reduce to the already known values in the various limits for 
the interpolating gauge parameters. In particular, as a useful check, it is verified that the beta 
function of the gauge coupling is independent of all the gauge parameters. In the concluding 
section we underline a potential application of the interpolating gauge for lattice simulations. 
Some technical details are collected in Appendix A. 

*In the case of the maximal Abelian gauge, the index a refers only to the N(N — 1) off-diagonal generators of 
SU(N), while the index A runs from 1 to iV 2 - 1. 
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2 The interpolating gauge fixing term. 



Let be the Lie algebra valued connection for the gauge group SU(N), whose generators T A , 
satisfying T A ,T B = if ABC T c , are chosen to be hermitean and to obey the orthonormality 

condition Tr (T A T B ^j = T F 5 AB , with A, B, C = 1, . . . , (iV 2 - l). As needed in order to dis- 
cuss the maximal Abelian gauge, we decompose the gauge field into off-diagonal and diagonal 
components, namely 

A^ = A A T A = A a ^T a + A\T\ (2.1) 

where the indices i, j, k label the N — 1 generators of the Cartan subalgebra. The remaining 
N(N — 1) off-diagonal generators will be labelled by the rest of the small Latin indices. For 
later use, we recall the Jacobi identity 

j ABC fCDE _|_ jADC jCEB + jAEC fCBD _ q ^2) 

from which it can be deduced that 

jabi jbjc + jabj fbci = Q) (2.3) 

jabc jbdi _|_ jabd jbic _|_ jabi jbcd g 

jabc jcde _|_ jabi jide _|_ jade jceb _|_ jadi jieb _|_ jaec jcbd _|_ jaei jibd q 

The field strength decomposes as 

F, u = F A T A = F* u T a + Fj lv T i , (2.4) 
with the off-diagonal and diagonal parts given respectively by 

F% = DfAl-DfA^+gf^Al, (2.5) 



F%, = d^Ai-d v A^ + gf M A^Ai, 



where the covariant derivative is defined with respect to the diagonal components A^ 

D ab = d ^ab_ gf aU A i^ (2f}) 

For the Yang-Mills action one obtains 

Sym = -\J d A x (F- u F^ a + F; U F^) . (2.7) 

In order to write down the interpolating gauge, we shall consider the most general gauge fixing 
term compatible with the decomposition of the gauge field into diagonal and off-diagonal com- 
ponents, according to eq. (|2.1j) . Let us first establish the nilpotent BRST transformations of the 
fields. From 

sA A = -D AB c B = -(d,c A + gf ABC A B c c ) , 

SC A = 9_ f ABC c B c C 
2 

sc A = b A , 

sb A = 0, (2.8) 
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one easily gets 



sAl = -(pfc b + gf abc Ay + gf aU A h /), sA^ = - (d,J + gf tab A«c l 



sc 



sb a 



gf abi c b c i + ^f abc c b c c 
0, 



sc 
sc' 



dfiabab 

2 



sb l = . 



2.i 



Here (c a , c l ) are the off-diagonal and the diagonal components of the Faddeev-Popov ghost field, 
while (c a ,b a ) and (c*,6 l ) are the off-diagonal and diagonal antighosts and Lagrange multipliers, 
respectively. We also observe that the BRST transformations (|2,9|) have been obtained by their 
standard form upon projection on the off-diagonal and diagonal components of the fields. 

Thus, for the most general BRST invariant gauge fixing term we write 



5gf 



(fx 



' 0M£ + — b a + a 2 gf abi A b A i n + c 4 ( + —b l 



+ ^gf ahc c a c*(f + ^-gf abi c fl c b c l + ^-gf abi c a c b c i 



(2.10) 



where cti, a 2 , «3, «4, as, a$ are six independent parameters. It is not difficult to see that 
the gauge fixing (|2.10|) interpolates between several gauges, for instance: the linear covariant 
gauges, including the Landau gauge as a special case, the maximal Abelian gauge, the modified 
maximal Abelian gauge and the Curci-Ferrari gauge. This will be checked in the next section, 
where the various limits for the gauge parameters will be analyzed in detail. 

The expression (|2.1()j) is easily worked out and yields 



Ql 



b a [ + ^b a + a 2 gf aU A\A^ + ^gf abc c b c c + ^gf abl tc 



«4 f abc-b r 



^5 rabi—b A 



(fx 



+ -fgfc^j + b l [d^A^ + -j-b 1 + -jdf c a c ) + c^V + cW 
+ gf abc d»(?c b A^ + g f - gf^d^dA^ - gf abi c a d> M c i A b l 



+ a 2 gf abt c <l d IM c b A l IM + a 2 gf abl (?d»c i A b 
+ a 2 g 2 f abi f bcd c <l c d A^A i > M + a 2 g 2 f abi f bcj (? d A c ^ + a 2 g 2 f abi f cdi (fc d A c ll A b ^ 

+ ( ?±g 2 f abc f adi c*<fc d C i + ^±g 2 f abc f ade ^c C C d C e + ^-g 2 f abi f cdi ^C C C d 



a 2 g 2 f abi f bcj c a c c A i u A jfl 



+ ^V^H^cW' + ^lg 2 fbijacd-b-i 



c c c d 



(2.11) 



For calculational purposes it is useful to eliminate the Lagrange multipliers. This is done by 
means of the equations of motion 



8 (Sym + Sgf) 
5b a 

6 (Sym + Sgf) 
5ti 



eabi tb 



«4 nta bc-b r 



a 5 „ rabi-b J , a 6 n rabi b - 



ai b a + dM£ + a 2 gf abl A b u A^ + ^gf abc -(fc c + ^-gf abl c b c l + ^-gf abt c b c 



azV + d' 1 A i fl + ^-gf abi c a c b . 



(2.12) 
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Substituting the expressions (j2.12j) into the expression (J2.1H) . we find 
Sgf = J d A x 



— —dPAIVAZ - —d^Aid u Ai - —gf abi d^AlA b y A iu 



2 

- ^g 2 f abi f acj A b l A i ^A c v A^ + c^V + gW + —gf abc c a d fl c b A'; l 
2ai ^ " 2ai 



+ 



+ 



2ai 

a 

V2ai 



f — + l) gf abc d^c b Al + (g- + l) gf abi d»<fc 



<?c b A^ + f ^ + l) gf M Ay^ 



— + gf^d^jAl - ^gf abi d^c a eA b + ( — + gf abi c a d fl c b A i , 
2ai / 2ai V2a 3 / * 



+ «2 



— ) gf M ?& i t*A b lt + a 25 2 ff abi f bcd + ^-f«bdfbci\ ^ c d A c^ 

— + l) a 2 g 2 f abi f bcj ^d'AlA^ - ^^g 2 f abi f ac i c b ? A C ^A^ 
a25 2 / aM / fecj c a c c ^yl J > + a 2 g 2 f abi f cdi ^c d A c ^A b ^ + ^ 5 2 



a 2 ( a 2 
"3" "5/ J + "4/ J H J J 

8 \ a 3 



+ ^ fabdjacA ^c c d c i + 

ai / 

+ ^Af^fde j ^cVc^ + Y + g 2 f abi f acj ^^c c c i + ^-g 2 (j M f 



^±jabcjadi\ ^ 5 i c c c d 

ai / 



(2.13) 



3 Interpolating among various gauges. 



Let us now proceed by showing how expression 1)2. 13|) can be used to interpolate among various 
gauges. We begin with the linear covariant gauges, including the Landau gauge. 



3.1 Linear covariant and Landau gauge. 

The linear covariant gauges (LCG) correspond to the following gauge fixing term 

S LCG = sJd 4 xc A (^b A + d^A^ 



which yields 



c A d»D AB c B + b A 



Slog = j d A x 
Eliminating the Lagrange multiplier 

S LCG = J (fx (c A d»D AB c B - ^-d^A A d u A A 

As usual, the limit a — > 0, defines the Landau gauge. In order to recover the linear covariant 
gauges from expression (j2.13|) . the values of the gauge parameters should be taken as is shown 
in Table 1, namely 

02 = CH4 = «5 = «6 = , (3-4) 



(3.1) 
(3.2) 

(3.3) 
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gauge parameter 


a\ 


02 


a 3 


a.4 


a 5 


«6 


limit 


a 





a 


a 6 


a 6 






Table 1: Values of the gauge parameters for the linear covariant gauges. 



and 



After this, expression (|2,13|) yields 



a\ = a 3 = a 



S*GF 



1 



1 



cPx ( -— d»A°;,d v Al - —d^Ai,d u Ai + c a d 2 c a + cW 



2a" ' » v 2a * 

•abcf^u,-ab ac , „ rabi au--a Jb \i i „ fabi Aab^n-i 



(3.5) 



"'"';>/'„"„' I''' \ ^3 g) 



+ gf aoc d fl & 1 d'A c fl + gf Mn d^c°A^ + gj am A^d^c 1 - gf^d^c'A^ 

Now, collecting together both diagonal and off-diagonal indices, we can rewrite the above ex- 
pression as 

5 GF = J d'x (-±d^A A d^A A + c A d 2 c A + gf ABC d»c A c B A c ^ . (3.7) 

Upon integration by parts of the last term, we obtain 

5 GF = J d'x {-±-d»A A d»A A + c A d»D AB c B ^j , (3.8) 

which is precisely expression (|3.3|) . In particular, the Landau gauge is thus recovered in the limit 
a -> 0. 



3.2 Curci-Ferrari gauge. 

Let us consider now the Curci-Ferrari gauge. It is a nonlinear covariant gauge, introduced in 
HUGHE], specified by 



S CF = s Jd 4 x (c A (^b A + d»A; 



-9f ABC c B c C 



(3.9) 



namely 



S CF = fd A x [b A ( -b A + d»A A - - g f A ^c"c u ) + c A d»D AB c* - -g A f A " u f uu *c A c»c u c* ) , 



o 



ABC^BC \ , -AwnAB„B a 2 f ABC fCDE-A-B n R 



X " 2 

which, after elimination of the Lagrange multiplier, becomes 



(3.10) 



S C f = f d A x(-^A A d v A A + c A d»D AB c B + 9 -gf ABC d»A A c B c c 



a 



(3.11) 



where use has been made of the Jacobi identity (|2.2j) . To obtain the Curci-Ferrari gauge from 
expression (|2.13|) the gauge parameters have to be taken as in Table 2, i.e. 



a 2 = , 

C*l = 03 = —04 = — CK5 = — «6 = a 



(3.12) 
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gauge parameter 




«2 


a 3 


«4 


a 5 


a 6 


limit 


-Q 6 





-"6 


«6 


a 6 


— a 



Table 2: Values of the gauge parameters for the Curci-Ferrari gauge. 



This leads to 
•Scf 



= / d A x (-—d»A a n d v A a v - —d^Aid u Al + (fd 2 c a + cW - -gf abc c a d^c b A c u 
+ ^f^d^Al + \ 9 r U d^c h A^ + ±gf M Al<*&><* - -g^d^A^ 



_ ^2 / yabc yacii j?abd y ac A C^C^C^C? OL-— ( j ^ _|_ y a ' )e ^cde jaci jbdi 

4 V / 8 V 

_ jacejbde^ ga^c^ _ jafei jacj-fe-^c^ _ «^2 ^jabi jacd + jabcjadi^ . 

(3.13) 

Collecting the various indices and using the Jacobi identities (|2.3JI for the quartic ghost interac- 
tion terms, we obtain 

S CF = f d A x (-—d»A A d v A A + c A 8 2 c A - -gf ABC c A d»c B AC + -gf ABC 'd"c A c B A° 
J \ 2a 2 yJ 2 yJ * 



^ „2 rabi facd-jb-c d J, , „ 9 ( rade febc i radi fibc\ -=a^b„c„d ® 2 rabi racj -Jb-i „c A 

■ c c c c — -g j j J c c c <y 



_ a g 2fahdjaci^ c c c d'} 



(3.14) 



In other words 



Scf = J d A x (-±-d»A A d»A A + c A d\ A - \gf ABC c A d^c B A^ + ^ f EAD f 



2 f EAD f EBC-A-B c C c D 



+ ^f ABC d^c A c B A^ 



(3.15) 



(3.16) 



Finally, integrating the last term of the above expression by parts, eq. (|3,15j) . yields 
Scf = j d*x(-^A A d v Ai + c A d»D AB c B + ^gf ABC c A c B d»A° 

+ ® g 2 fEAD f EEC d A-B c C c D 

which reproduces the Curci-Ferrari gauge fixing term, ()3.11|) . 
3.3 Maximal Abelian gauge. 

The so called MAG gauge condition amounts to fixing the value of the covariant derivative, 
D^A^, of the off-diagonal components by requiring that the functional 



n[A] = (vt)- 1 J d 4 x [a^a^ , 



(3.17) 



attains a minimum with respect to the local gauge transformations. This corresponds to impos- 
ing 

However, this condition being non-linear, a quartic ghost self-interaction term is required for 
renormalizability purposes. The corresponding gauge fixing term turns out to be [9*1 110] 



(3.19) 



and the real MAG, ea. Q3.18j) . is obtained by setting a = after renormalization. The MAG 
condition allows for a residual local ^(l)^ -1 invariance with respect to the diagonal subgroup 
of SU(N). In order to have a complete quantization of the theory, one has to fix this Abelian 
gauge freedom by means of a suitable further gauge condition on the diagonal components A^ 
of the gauge field. A common choice for the Abelian gauge fixing, also adopted in the lattice 
papers [l"QEl, is the Landau gauge, given by 



S" 



diag 
MAG 



s I cPx&dpA* . 



The complete gauge fixing is, thus, the sum of expressions (j3.19]l and (|3.20j) . namely 



Smag = sjd 4 x Y [DfA h » + ^b a j + c%A^ - ^gf^c'd - ^ g f abc c a c b c c 
Working out expression (]3.21j) we find 



(3.20) 



(3.21) 



Smag 



OL , 



a 



b a (d^A* - gf^A^A 1 " + ^b a - agf abl (fc l - ^gf abc c b c c j + 6*0" AJ, 

+ <?D ab D bc »c c + cW + gf^A^d^c 1 + g f abi A^ d» j + gf abi A a ^c b d^c l 
+ gf^d^cf^A^ - g 2 f abi f bcj c fl c i A^A j ^ - g 2 f abi f bcd <f c d A£A^ 
- g 2 f abi f cdi c a c d A b A Cf * - -g 2 f abi f cdi c a c*c c c d - -#*f abc f adi &Z s d i j 



-g 2 f abc f ade c*c c c d c e 



(3.22) 



The task of eliminating the off-diagonal Lagrange multipliers is straightforward. To eliminate 
the diagonal Lagrange multiplier we have to introduce a parameter (3 via the term ^6*6*, and 
then let (3 — > to recover the transversality of the diagonal gauge field. This will lead to the 
following gauge fixing 



Smag 



—d» A a u d v Al - —d^A„d u Ai + —f ahi d^A^.A^,A i 
2a 2/3 a " v 



lim d 4 x 

a 2 I 

!L_fabifacj A b^ A c A jv + <fi D ab D bc» c c + Jtflj + g jabi jaj Q^i + jabcQp^Jj^ 

^ g jabc-gi Qfi J} A c 9_ f<2 jo-bi jbcd jabd jbci\ A c jji^i ^2 jabi jcdi^a^d pb A c^i 



fabc fadi\ ^^c^d^i 



^Lg^ ^2 f bc * jdei _|_ jabc jade jabd jace^j (^(^(^(f 



(3.23) 
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gauge parameter 


a.\ 


a? 


a 3 


CtA 


a 5 


a 6 


limit 


a 


-1 


/3 end () 


—a 


-2a 






Table 3: Values of the gauge parameters for the maximal Abelian gauge. 



Expression (|3.23|) is recovered from the interpolating gauge fixing (|2.13|) . once the gauge param- 
eters are taken as in Table 3, 



"2 



and then 



«5 

— ola = = a 

2 

- 1 , a 6 = , 
a 3 = p _> . 



(3.24) 



(3.25) 



Therefore, expression (|2.13j) becomes 



lim I d 4 x 

0^0 



—d»A"d v A a v - —d^AiWAl + -gf abi d»AlA b v A l 
2a * v 2(3 a " ' 



1 



1 

2' 



+ ^gf abc d^c b A^ + gf abi d^c b Ai + gf abi Alc b d^ - g f abi d» c b A^ 

2 | jabi yftca! _ \ jabd J bci \ c a C d A C A ip " + g 2 / afc * f bc ^C a C C A i A^ — g 2 f abi f cdi C a C d A C A b ^ 



a 



^2 / yafec jadi ^ jo-bd J ac ^ (fc C C d C^ 



q 9 f 2 jabi jcdi _|_ jabe jcde jace jbde^ (^(f q c q" 



(3.26) 



which, after rearranging some indices, yields the expression (|3.23|) . 
3.4 Modified maximal Abelian gauge. 

For completeness, let us include in the present analysis also a slightly modified version of the 
maximal Abelian gauge. This gauge, discussed in |13j . is obtained by making the following 
choice for the diagonal gauge fixing term 



SmIg = sJd A x (c^4 + gf ab ^A b ^) , 



(3.27) 



Thus, for the complete gauge gauge fixing in the modified maximal Abelian gauge, MMAG, we 
have 



SMMAG = 8 J d 4 X (V (&>A a u + % a ) +7fV^ - ^gf^cfj - ^gP bc c a C h t 



(3.28) 



namely 
Smmag 



/ 



d 4 x 



Ot , 



b a [ &*A^ + ^b a - agf^cfj - ^gf abc c*c c ) + + (?d»D ab c b 
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gauge parameter 


a\ 


«2 


a 3 


«4 


a 5 


a 6 


limit 


a 





f3 end Q 


—a 


-2a 






Table 4: Values of the gauge parameters for the modified maximal Abelian gauge. 

+ cW + gf^A^d^c 1 - gf abc d»(?Alc c - g^d^A^c* - ^-g 2 f bc f^^^d 



-g 2 f abi f cdi ^c b c c c d - -g 2 f abe f cde c a ^c c c d 



which, upon elimination of the Lagrange multipliers, becomes 



(3.29) 



Smmag = lim / d A x 



^-d^A a u d u A a u - ^-d„A ilx d v A iv + <?d»D ab c b + cW 
2a F 2p ^ 



+ gf aM A a ^c b d^ - ^f abc ^d^c b Al + ^f abc d'"c a c b A c fl + gf^A^j 



—g^ ^2 _|_ jabe jcde _|_ jade jbce\^ -a-b^c^a 



(3.30) 



Again, expression ()3.30|) follows from ()2.13j) by setting the gauge parameters as in Table 4, 

"5 



Q'2 



and 
Thus 



— C*4 : 

Z 

"6 = 0, 

a 3 = (3 _> . 



(3.31) 



1 

2 



1 



5 GF = hm / d 4 x ( -J-aM^S"^ - -L^A^Al + cfdV + gW - ±-gf abc (fd fl c b A c tl 



- ^gf abc d^c b Al + gf abi d^c a c b A i fl + gf aM ^c'W + gf abi ?d>*(?A b tt 



^ ^afcc yacii 2 jo-bd j ac ^j c b C C C d C^ 

2 

q 9 ^ johi jcdi _|_ yafee ycde jace jbde^j 



(3.32) 



By rearranging the indices of the quartic ghost terms by means of the Jacobi identities (|2,3|) . it 
is then easy to show that 

Sgf = 5mmag • (3.33) 



4 Multiplicative renormalizability of the interpolating gauge. 



To prove renormalizability, let us first write down the Slavnov- Taylor identity corresponding to 
the BRST invariance of the classical action. Following the algebraic renormalization setup |14j . 
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we introduce a set of external sources coupled to the nonlinear BRST transformations (j2.9|) of 
the fields, 

5 ext = sjd 4 x (-VL^A^ - fi^Aj, + L a c a + LV) , 



with 



= sCT u = sL a = sL J = 



(4.1) 
(4.2) 



or 



S cxt = J d^x [-W (Dfc b + gf abc Alc c + gf M A b ^) - fi* (fyt? + <7/ o6i A£c< 



rabijbji i 5 fabcb 



+ L a gf am c b c l + ^f abc c b c c + ^f abl c a c b L 



9 tabi „a„b i 



Therefore, it follows that the complete action X 

£ = 5Vm + Sgf + Sext ; 

obeys the Slavnov- Taylor identity 

<5E 5£ <5£ «JE «5E SE 6E £E 



5(E) 



+ 



H 1 r r + f 1 h 6* r 

6L a <5c a 5L l 5c l 5c a 5c 1 



. 



(4.3) 



(4.4) 



(4.5) 



In order to characterize the most general invariant counterterm which can be freely added to any 
order of perturbation theory, we perturb the classical action £ by adding an arbitrary integrated 
local polynomial E co in the fields and external sources of dimension bounded by four and with 
zero ghost number, and we require that the perturbed action (E + ^E count ) satisfies the same 
Ward identities as E to the first order in the perturbation parameter rj, being, 



5(£ + 7/E count ) = +0(t? 2 ) 

This amounts to imposing the following conditions on E count 

5sE count = q ^ 

where 5x is the linearized Slavnov- Taylor nilpotent operator, 



(4.6) 



(4.7) 



5 E 



<JE 5 SE 5 <5E 5 



<5S> a 5Af, 5 A" 5Q^ a Stt^ 5 A], 



+ 



<5£ 5 (JE i 
+ 



8A\ 5L a 5c u 



<5E 5 <5E 5 5E 5 6 ri 5 
+ 1 r — H r — r + b a Yb r 

5c a 5L a 5L % be 1 5c l 5L l Sc a 5c l 



(4.8) 
(4.9) 



5 S 5 S = . 

From the conditions (|4.7j) . it turns out that the most general invariant counterterm E count can 
be written as 

E count = a S YM + SxA- 1 , (4.10) 
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A 


c 


c 


b 


n 


L 


dimension 


1 





2 


2 


3 


4 


ghost number 





1 


-1 





-1 


-2 



Table 5: Dimensions and ghost numbers 



where A 1 is a local polynomial in the fields and sources, with dimension four and ghost number, 
which reads 



/ 



(fx I ai ^b a (? + a 2 a 2 gf abi <?A b u A i » + a 3 ^6V + a^gf abc c a c h c c 



+ a^gf^c^d + a^gf abi c a c h c l + a 7 ^A»c a + A 1 / + a 9 ^^ + a 10 %A^ 
+ anL a c a + a i2 LV) , (4.11) 

where ai \ i € {0, 1, . . . , 12} are independent parameters. These parameters can be reabsorbed 
by means of a multiplicative renormalization of the gauge coupling constant, of the gauge pa- 
rameters, and of the fields and sources, according to 

S($ o ,0o,«/o,io,eo) = S(cI>,0,J,i,O+^ count (cI», ( /),J,i,e) , (4.12) 

where 

$o = Zl /2 <Z>, 

4>o = z l J 2( t> , 

Jo = ZjJ , 
jo = Zjj , 

£o = (4.13) 



with 



$ = {A a ,b a ,c a ,c a } , 
<P = {A\b\c\?} , 

j = {n a ,L a } , 



j 



£ = {g,a 1 ,a 2 ,a 3 ,a4,a 5 ,a 6 } . (4.14) 

For the renormalization of the off-diagonal and diagonal gluons, and of the coupling constant g, 
we find, respectively, 



7-1/2 


1 _L 


7I/2 


1 _L ( ^ 




1 ao 



(4.15) 



12 



while the Lagrange multipliers renormalize as follows 

ZJ = I + 77 \a 7 - — 



2 



Zl' 2 = l + v(a 8 - a f). (4.16) 

The renormalization factors of the gauge parameters are 

Z ai = 1 + r] (ct + ai - 2a 7 ) , 

Z a2 = 1 + r] (a 2 - a 7 ) , 

Z a3 = 1 + 77 (a + a 3 - 2a 8 ) , 

Z a4 = 1 + r] (a + a 4 - 2a 7 ) , 

Z a5 = 1 + 77 (a + 05 - 2a 7 ) , 

Z afi = 1 + 77 (a + a 6 - a 7 - a 8 ) . (4.17) 

The renormalization of the ghost and anti-ghost fields are constrained by the following self- 
consistent relations 



zl l2 zl' 2 = 


1 + 77(07 


- an) 


zl /2 z l J 2 = 


1 + 77 (a 8 


- au) 


zl /2 zl' 2 = 


1 + 77(07 


- a 12 ) 


zy 2 zi' 2 = 


1 + t? (a 8 


- a u ) 



which give 



(4.18) 



rri/2 1 , ( a 7 a 8 an a 12 

= 1 + _f£7 08 aii_ai2 
+ ' V 2 2 + 2 2 
7 i/2 /a 7 a 8 an 012 

Zc " " 1 + Hy~T ~~~T~ 

4 1/2 = i + ^(-y + y-f -f ) • (4-19) 



Finally, for the renormalization of the external sources 

^ -, (a 7 a 8 a u a 12 \ 

Zn = l-v[ Y + Y + a 9 + — + — j 

^ -, f a 7 a 8 an . a 12 

Z L = 1 + 7/ - a 7 - a 8 - ai2 

= l + 77(y-a 7 -a 8 -an) • (4.20) 
This completes the proof of the multiplicative renormalization of the interpolating gauge fixing. 
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5 One loop renormalization. 



We now record the details of the explicit renormalization of the interpolating gauge at one loop 
which requires the computation of the renormalization constants Za, Za-> %ibi Zgi Z c , Z z , Z c , Z g 
and Z ai , with i = 1, . . ., 6, and where Z^ is the quark wave function renormalization constant^. 
Given the explicit form of the Lagrangian we do not need to consider the renormalization of all 
three- and four-point interactions. The procedure we adopt is to determine the renormalization 
constants of the fields and of the gauge parameters 01,03 before deducing the explicit values of 
the remaining four parameters 02,04,05,06- We use dimensional regularization in d = 4 — 2e 
and subtract the divergences using the MS scheme. In previous work we carried out similar 
renormalization in the Curci- Ferrari gauge and MAG, |15l 1161 HI I17j. using the Mincer algo- 
rithm, |181 119j . written in Form, [2U], where the diagrams for the appropriate Green's function 
were generated by the Qgraf package, |2J. However, to determine the full interpolating gauge 
renormalization constants it transpired that we needed to renormalize several four-point inter- 
actions in addition to various three-point functions. In this instance the Mincer algorithm is 
not fully appropriate since it can only be applied to two-point functions with massless propaga- 
tors. For a four-point function one would therefore need to nullify two external momenta which 
could potentially introduce spurious infrared divergences. To circumvent this difficulty for the 
four-point function renormalization, we applied an alternative algorithm, still in dimensional 
regularization, which involved nullifying all external momenta but systematically introducing 
a temporary infrared mass regularization according to the approach of |221 I23| . This involved 
repeatedly replacing massless propagators by massive ones via 

1 = 1 + (2kp - p 2 - m 2 ) 
(k — p) 2 (k 2 — m 2 ) {k — p) 2 (k 2 — m 2 ) ' 

until the contributions from the last term are finite and can be neglected by Weinberg's theorem, 
|24| . The resulting one loop vacuum bubble Feynman diagrams are then elementary to evaluate. 
For the three-point functions we will still nullify an external momentum but we will do so 
only in Green's functions where we know no infrared singularities will be introduced. In either 
approach we will still use the method of to extract the explicit renormalization constants. 
This involves computing the Green's functions in terms of bare quantities before transforming 
to renormalized variables and parameters. The form of the tree term then allows one to fix the 
renormalization constants associated with that Green's function. This approach avoids the need 
to consider subtractions and moreover is more appropriate to performing automatic symbolic 
manipulation computations. One additional point arising from using the method of [2^j is that 
we cannot renormalize 05 itself since it never appears as an isolated coefficient in one type of 
interaction. However, since it appears in the combination (1 + 05/(201)) we have chosen to 
define a new parameter 

a 5 = 1 + P~ , (5.2) 
2oi 

and determined its renormalization constant Z„ 5 . Clearly this choice does not invalidate the 
renormalizability of the interpolating gauge. We note that to recover the various gauges discussed 
in Sect. 2 the value of 05 becomes 1,^,0 and respectively in the linear covariant, Curci- Ferrari, 
maximal Abelian and modified maximal Abelian gauges. 

^ Although we did not consider matter fields in the previous sections, it can be easily checked that the inter- 
polating gauge fixing term 12.1311 remains renormalizable if spinor fields are included. 
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Since we are considering a gauge where the color group is split into diagonal and off-diagonal 
sectors we briefly recall the essential properties of the group algebra which is discussed in more 
depth in [T7j . In addition to the usual color group Casimirs C F , Ca and T F for a general color 
group we define the dimensions of the centre of the group as iV^ and the dimension of the 
off-diagonal sector as iV|. Clearly 



n£ + NX 



N A , 



(5.3) 



where Na is the dimension of the adjoint representation. For instance, for SU(N) one has 

that 

NX - 2Ni C A , f icd f Jcd = C A 6 ij , 



NX = N - 1 and NX = N(N - 1). Consequently, [Hj, one can show that 



eiab fiab 



jiab ji 
jacj jbcj 

japq jbpr jcqr 

jipq jbpr jcqr 



rry I rparpb 



nXCa , f abc f abc 
N d 



'-^-C A S ab . 

NX 

[NX ~ 3JVj] 
2NX 

[NX-^Nj] 
2N° 



jacd jbcd 



c A f 



abc 



apq fbpi rcqi 



japq jOpi j 



N d 
A /~i rabc 

2N£ AJ 



Cj^f ibc f ipq f b Pif c Qi = —A Q A f 



ibc 



rpbrparpb 
rj-iirpCLrpi 
rparpirpa 



T F 6 ab , 




= 


1^ 


, T a T 


a, 


— 


C F - 


C A 
2 


—NX 

Np A 


—NX 

N F A 


C A Nf 
2NX J 




—NX 

N F a 


C A ' 
2 _ 


rj-ii 


5 



Tr (VT^ 

c a n£ 



o 



T F 5 ij 



2NX 



rpa 

1 1 



rpjrplrpj 



Tf 
N F 



N d T t 



where N F is the dimension of the fundamental representation and satisfies 



C F N F 



m + m 



T F . 



(5.4) 



(5.5) 



For completeness we note that we have assumed that f acd f bcd \ s proportional to 6 ab which is 
certainly true in SU(2) and we have checked that it is valid in SU(3), Moreover, we have 
assumed that T l T l is proportional to the unit matrix, /, which is certainly true for SU(N). 

Before turning to the explicit discussion of how we renormalized the interpolating gauge 
Lagrangian we define the anomalous dimensions of the wave functions and parameters in terms 
of the renormalization constants we will compute. We have 



d_ 
da 



) d 
da n 



7ai W 



1 — a\ 



-k— In Z a 
oai 



1 -l 



d 

7 A (a) + 0(cl)q^ + 



E 



a«7a„ (a) 



d 



InZ, 
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Ja 2 (a) 

7a 3 («) 
7Q4( a ) 

7s 5 (a) 
7« 6 («) 



1 + CK2 




1 - a 3 



a — In Z az 
da* 



-1 



7a(o) + /3(a) 



<9a n 



lnZ 



d_ 

da 



1 + 04 



1 + a 5 



1 + a 6 



d 

— — In Z aA 

00!4 



1 -1 



/3(a) 



9a 



+ Xl a ™7an( 



n' 

8 



Q-2 ) 



) — 

3a r 



lnZ, 



03 > 



-1 



-1 







) — 

da n 

d 

da n 



lnZ, 



«4 ) 



lnZ, 



Q5 ' 



lnZ n 



(5.6) 



where Yn' nieans a sum over 1 to 6 but where the term corresponding to the parameter 

associated with the left hand side is omitted. We have defined a = g 2 /(16tt 2 ). Moreover, 

4> € {^4, A 1 , cc, (?c l , <?c, cc l , tp} where we now use the notation that a superscript i on a field 

indicates that it is a diagonal object. We have chosen to define the anomalous dimensions of the 

ghost and anti-ghost renormalization constants in both sectors in a non-standard way. This is 

because in the general interpolating gauge to extract individual anomalous dimensions for each 

II 

field itself entails making an extra assumption. So, for example, Z 5c = Z? Zc ■ We also note 
that our convention for defining the relation between the bare and renormalized parameters for 
cti and «3 is 

"io = Z~lZ A a.\ , a 3o = Z~jZ A a 3 , (5.7) 

where denotes the bare quantity. 

We now present the actual details of our renormalization. In order to extract all the renor- 
malization constants we have had to select a particular set of Green's functions to renormalize. 
Since the order in which we consider the particular interactions is crucial, it is worth refering 
to the actual Feynman rules given in the Appendix A. First, applying the Mincer algorithm, 
|18|I19|. to all the two-point functions we have already discussed, it is straightforward to extract 
the wave function renormalization constants and those for a\ and a 3 and encode them in the 
anomalous dimensions as 



7a(o) 



1 



6N; 



-[NX(C A (3 ai -13) + 8T F N f ) 

+ N£(C A (-3ai - 9a 2 + 3a 3 - 3a 2 a 3 ))l a + 0(a 2 ) , 



j A t (a) 



[C A (3aia 2 + 3ai + 9a 2 - 13) + 8T F N f ] a + <3(a 2 



lai ^ = \2aiN° N ^ ( C '^( _6a i ~~ 6 «i Q 4 + 26ai - 3a 2 ) - WaxTpN^j 



+ N A C A (l8a 2 a 2 - 12a?a 5 + 18a 2 - 6aia 2 a 3 - 18aia| 
— 12aia 2 a 3 + 6ai« 2 a6 + 18aia 2 — 6aia 3 

— 2 2 

+ Yla\a^ — 12a\a^ae + 6«ia6 — 36a 2 + 6a 4 
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7a 3 (°) 
7cc(a) 



+ 0(a 
1 

1 



Ca (— 6aia2«3 — 6ai«3 — 6a2a3«6 — 18a203 — 60306 + 2603 — 3a, 



16a 3 T F N f ]a + 0(c 



AN; 



7cv( a ) 
7^ (a) 



- CU(a a - 3) 

+ iV^C74(3a:ia2 — 2a\a^ + a\ + a\a 3 — 3a\ — 2a2«3 
+ a 2 «6 + 3«2 + «3 + "6)] a + 0(a 2 ) , 

«i — Q2Ct6 — 3«2 — «6 — 3] a + 0(a 2 ) 
[NXT F ( ai - a 3 ) + N F C F a 3 ] a + 0(a 2 ) , 



[— ot\Oi2 + 2aia5 

1 



(5.8) 



where Nj is the number of quarks. Due to our particular definition of the renormalization of 
the gauge parameters a\ and 03 we note that for completeness their anomalous dimensions are 



7 A (a) + 7«i(a) 



C A 



j / o 9 9 9 

I 6a 1 «2 — 4a 1 a5 + 4a x — 2aia 2 a3 — Qa\a 2 — Qa\oi2a 3 
+ 2aia2a6 + 4ai«4 — 4aia5a6 + 2aiQ6 ~~ 12«2 + 204 



7a* (a) + 7a 3 ( a ) 



4«3 



+ iV| (-2aia 4 - a|J a + 0(a 
2a 2 a 3 + 2a 3 + a 6 ] a + 0(a 2 ) . 



(5.9) 



The next stage is to consider the vertex renormalization. Equipped with the diagonal and off- 
diagonal gluon and quark anomalous dimensions we have renormalized the triple off-diagonal 
gluon, quark off-diagonal gluon and quark diagional gluon three-point functions and verified 
that the correct one loop QCD /3-function of emerges, independent of all six gauge fixing 
parameters, as 
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-C A - -T F N f 



+ 0(a 3 



(5.10) 



This is a useful check on our procedure since it is known that the /3-function is gauge independent 
in mass independent renormalization schemes. For the four remaining renormalization constants 
those for 02 and 04 are the easier to extract. From an examination of the Feynman rules then 
considering the respective vertices (A l A l cc) and {Acc} each parameter appears as a factor and 
therefore their renormalization can be extracted by the method of [25] given that Z 5c has been 
determined from the off-diagonal ghost two-point function and Z g has first been established 
as the usual value. In this instance for the three-point function the momentum flow has been 
carefully chosen so that there is no external momentum flowing into the c leg. We find 



7« 2 (a) 



C A 



8a 2 N£ 



N% (2a>ia>2 + 2ai«2 + + 602 



+ m 



20102^5 + 8ai«2 — 2aiQf5 + 2ai 



6«2 + o| a 6 + 2a2«6 + 6«2 + «6j a, + 0(a 2 ) 



7a 4 ( a ) 



1 



C A (3al + 26a 4 ) - 16a 4 7>iV / ) 
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+ N^Ca (l8aia2«4 — 12ai0206 + 12010506 — 12aia6 — 180^04 

o 9 \ 

+ 36o 2 + 602040:6 + 180204 — 12o 4 + 60406 J 



a + 0(a 2 ) 



(5.11) 

To extract Z a6 one first requires the renormalization constant Z 5 i c since the vertices where a$ 
appears always involves the combination of fields c 1 and c. To do this we consider the vertex 
(Affc) and choose the momentum flow through the off-diagonal gluon leg to be zero. From the 
Feynman rule this means that the tree term of the Green's function does not depend on any 
parameters Oj and therefore the only renormalization constant which is undetermined is Z 5 i c . 
Hence, 

7c» c (a) = [iV|C' J 4(-aia2 + «i — a2«6 - 3«2 + «4 — «6 — 3) 



+ i\^CU(-2ai + 2a 3 - 2a 4 + 2a 6 ) a + 0(a 2 ) . (5.12) 



Equipped with this we can then determine Z a(j in two ways. The first is to renormalize the 
same vertex again but nullify the momentum flow through the <? leg which leaves a tree term 
involving the as yet undetermined Z a6 . An alternative way is to consider the (AA l c l c) vertex 
given the diagonal and off-diagonal gluon wave function renormalizations and Z a2 have already 
been determined. We have considered both and deduced the same renormalization constant for 
Z a(i for both Green's functions which is a useful check. We found 

7a 6 («) = 12 ^ ro i N A (C A (-3aia 2 - 3a 2 06 - 9a 2 + 6a 4 - 3a 6 + 26) - 16T F N f ) 

+ N^Ca (9a\a2 + 60105 — 9oi + 302O3 — 9o2 + 60203 + 902Q6 

+ 9a 2 + 3a 3 - 12a 4 + 9a 6 )] a + 0(a 2 ) . (5.13) 

Finally, we follow a similar procedure to extract Za 5 - However, given we have just found Z ae 
by considering the vertex (ccc^c 1 ) then Zq, 5 can be extracting immediately given that the wave 
function renormalization for this vertex can be written into the already determined combination 
Z 5c and Z 5 z c i. This gives 

7a B (a) = [NX (2aja 5 - 2a\ + af) 

+ N% (-4aia 2 a6 + 4oio 5 a 6 - 4oia 6 + 12o| - 2af) a + 0(a 2 ) . (5.14) 

However, to obtain the complete set of renormalization constants we still need to deduce Z 5c i. 
Again this can be extracted from two vertices which are (Ace 1 ) and (AA l cc l ) where there is no 
momentum flowing through the c % leg in the former case. We have renormalized both Green's 
functions to ensure that the same renormalization constant emerges and found in both cases 

IcAa) = -^-[^^(20105-01-04-3) 



+ A^4Ca(3oiQ2 — 2oi05 + 3oi + Q2O3 — 3Q2 — 20203 

+ o 2 a 6 + 3a 2 - a 3 + 2a 4 - o 6 )] a + 0(a 2 ) , (5.15) 
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which is a useful internal check. 

Having produced the anomalous dimensions for the general interpolating gauge we are now 
in a position to take the various limits of the four special cases discussed above in order to 
check with the corresponding known anomalous dimensions for these gauges. In comparing with 
known anomalous dimensions any agreement we indicate is for all the anomalous dimensions 
except for the interpolating parameters 02, 0:4, 05 and a§. For the linear covariant gauge case 
we have 



7A(o) 
Ten (a) 
7 Q2 (a) 

7a 4 ( a ) 

7a 5 ( a ) 
Icda) 

7v( a ) 



lAi (a) = - (C A (3a - 13) + 8T F iV» a + 0(a 2 
6 



7a 3 ( a ) 
C A 



- - (C A {3a - 13) + 8T F N f ) a + 0(a 2 ) , 
N£ (a + 3) + 3Ni (a + 1)1 a + 0(a 2 ) , 



7Q6 (a) = - (13C A - STpNf) a + 0(a 2 
0(a 2 ) , 

7 glcl (a) = 7 5c »(a) = 7c* c (a) 
aCpa + 0(a 2 ) , 



(a - 3)C A 



a + 0(a 2 



(5.16) 



where for 7 Q4 (a) we have taken the limit 02 — > before setting 0:4 = 0. These expressions are 
in agreement with the well established linear covariant gauge anomalous dimensions. For the 
Curci-Ferrari gauge we have 



1a{o) 

7«i(o) 

7« 2 ( a ) 
7a 5 (a) 
7cc(a) 
7v(°) 



lAa) = g (CU(3a- 13) + 8T F iV / )a + 0(( 



1 



7a 3 («) = 7a 4 (a) = 7a 6 (a) = - ^ ( C ^( 3a " 26 ) + ™T F N f ) a + 0(a 2 ) , 



C A r 
8A2 



A^ (2a + 6) + A^ (5a + 6)1 a + 0(a 2 ) 

[a - 3)C A 



0(a 2 ) , 

7c*c*(a) = 7ec i (a) = 7c» c ( a ) 
aCpa + 0(a 2 ) , 



a + 0(a 2 



(5.17) 



which are also in agreement with [£l I27| I28| I29j . The equality of the four ghost wave function 
anomalous dimensions in this and the previous case is consistent with each of the individual 
wave function ghost anomalous dimensions being equal and equivalent to the value given. For 
the first of the MAG cases taking the limits with the values in Table 3 we find 



7a(o) 
7a* (a) 
7«! (a) 



1 



6N£ 



1 



3 



N£(C A (3a-13) + 8T F N f ) + NiC A (-3a + 9) a + 0(a 2 ) 
(HCa — 4T F Nf) a + 0(a 2 ) , 



12aN£ 



AT 



4 (cU(-3a 2 + 26a) - 16aT F N f } 
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+ NX Ca(-6o 2 - 36a - 36) a + 0(a 2 ) , 



la 2 {a) 

7a 3 (°) 
7« 5 («) 

lcc{a) 

7? c («) 
7^(o) 



0(a 2 ) , 



- {UC A - ATpNf) a + 0(c 



C A 



AaN 



12NX 



- (n£ (l2 - 2a 2 ) - a 2 iV|) a + 0(a 2 ) 



iYj , (C A (-3a + 35)-16rpA? f ) + A^C J 4(-6a-18) a + 0(a 2 



C A 



(fitf (-6 - 2a) +NX(a- 3)) a + 0{c 



0(a 2 ) , 
C A 



l -(Ni(-6-2a)-3NX)a + 0(a 2 ) , 



4 

aN£T F 



N F 



a + 0(a 2 ) , 
a + 0(a 2 ) . 



(5.18) 



These are in agreement with the one loop arbitrary color group anomalous dimensions recorded 
in Although unlike the previous two cases the combination of the ghost wave function 

anomalous dimensions are not equal, their values are consistent with the values given in |17j . 
Finally, we quote the values for the parameters stated in Table 4 are 



7A(a 
1 A i(a 
7ai(a 

la 5 (a 

7a 6 (a 
lcc(a 

7c C *(a 

7c* c (a 
7^( a 



6NZ 



NX {C A {3a - 13) + 8T F N f ) - 3aC A Ni a + 0(a 2 ) , 



1 



- ((3a - 13)C A + 8T F N f ) a + 0( 
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7a 4 ( a ) 
1 



12NX 



NX (C A (-3a + 26) - l6T F N f ) + 12aN£C A a + 0{a 2 ) , 



6 



((13 - 3a)CU - 8T F N f ) a + 0(c 



oQa 
AN 



±{2Ni + NX)a + 0(a 2 ), 



12NX 
C A 



NX {C A (-6a + 26) - l6T F N f ) + 3aN A d C A 



a + 0(a 2 



Am 



l -[aNi + NX(a-3))a + 0(a 



(a + 3)C A 



a + 0(a 2 



Ca 
AN 



l - (aNi - 3NX) a + 0(a 2 ) , 



A 

3C A 



a + 0(a 2 ) , 



aNXT F 
Np 



a + 0(a 2 ) , 



(5.19) 
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where 7 Q2 ( a ) is singular as ai — > 0. However, this singularity does not affect the anomalous 
dimensions of the gauge parameters entering the modified maximal Abelian gauge. We regard 
the results as new since as far as we are aware there are no similar computations in this 

specific gauge. 

6 Conclusion. 

We have constructed a six parameter family of covariant gauge fixings, which allows us to 
interpolate between various familiar gauges: the Landau, the linear covariant, the Curci-Ferrari 
and the (modified) maximal Abelian gauges. Using the algebraic renormalization formalism, we 
have been able to prove the renormalizability of the proposed gauge to any order of perturbation 
theory. In addition, we have also computed explicitly every anomalous dimension at one-loop 
order in the MS scheme, confirming that the obtained results reduce to the already known values 
in the various special cases, if already known. The gauge parameter independence of the gauge 
beta function was confirmed. 

We end by noticing that the concept of an interpolating gauge may also find potential 
application in lattice numerical simulations. More precisely, requiring that the functional 



with k the interpolating parameter, attains a minimum with respect to local gauge transfor- 
mations, leads to the MAG for k = and to the Landau gauge for k = 1. This could be 
useful in order to study the presence, in the MAG, of power corrections in to the strong cou- 
pling constant. Such power corrections are by now quite well established in the Landau gauge, 
[301 I3T1 OH ES]. As it turns out that one can interpolate from the Landau gauge to the MAG, 
one could therefore expect similar corrections in the MAG. One might also hope to obtain some 
information on the relevance of the Gribov copies in the MAG. 
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A Feynman rules. 

In this appendix we list the explicit Feynman rules we used in momentum space. For the 
propagators we have 





(A^p)A b v (-p)) 





PfiPu 
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(4(p)4(-p)> 

(c a (p)^(-p)) 



p 2 



P 



tab 



, (^(pJ^'C-p)) 



8W 

— , {i>{p)i>(-p)) 



pi, p. 

where p is the momentum. The non-zero 3- and 4-point vertices &re$ 

(4(pi)^(p 2 )V<(p3)} = gT%, 

(A a ( Pl )c b (p 2 )c c (p 3 )) = -igf abc 



p- 



(A.l) 



(A«(pi)c\p2)c k (P3)) = - ifjf 

(A«( Pl )ci(p 2 )c c (p 3 )) 

(Al( Pl )c b (p 2 )c c (p 3 )} 
(A^ Pl )Al(p 2 )A c a {p 3 )) 
(A;( Pl )Ai(p 2 )A k a (p 3 )) = igf 

{Al(pMl(P2)Al{p 3 )A d Jp A )) 



abk 



1 + ) Pi - P3 
2ai/ / /j, 



-"5Pi - (1 + a 2 )p 3 



^/ ab ° (^<r(P2 - P3)/x + ^(P3 - Pl)i> + V(Pl ~ P2)<r) , 



1 + 77~~ ) Pi - (1 - "2)P3 

2«3 / / M 



afefc 



P2 Pi - P3 + Va/j, [ P3 -\ P2 - Pi 

Oil J p \ OL\ /v 



+ Vpvipi -P2)u) , 

'Hpp'Hucr) 

+ fd dbC (VpvVcrp — V/MrVvp) + fo bCd {VpcrVvp ~ ^pp^va) 
Vp.p'Hvtj) + fo {VpisVcrp T lpcj T lvp) 



(A«(pi)A b u (p 2 )A*(P3)A l p (p4)) 
(A«(p 1 )4(p 2 )4(p 3 )4(P4)) 

(A;( Pl )A b u(P2)c c (p 3 )c d ( Pi )) 
(A;( Pl )Ai(p 2 )c c (p 3 )c d ( Pi )) 

(^(p 1 )4(p 2 )c C (p 3 )c'(P4)) 

(A a J P i)Ai(p 2 )c k (p 3 )c d ( Pi )) 



g 2 ^abcl (_ v ^ Vup + n^r,^) + f% cbl (-V^rjap + TlwVua) 
+ fo WC (-VpuVap + VpaVup) 

g 2 ( ff bl 



VpvVcrp ^ VpuVvp i VpcrVvp ~r 1}ppVv 



+ fo 



albk 



a 2 



g 2 (/r M «2^+/d ca %^ 

<? 2 (fZ l ^Vpu + f^a 2V ^ 



9' 



-Jo 2oi ^ 



'We note that in |17| the corresponding Feynman rules for the four-point corrections did not all have the 
correct coupling constant factor of g 2 and the {AAA 1 ) vertex was inadvertently omitted. 
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(4(^1)4^2)^(^3)^(^4)) = - 9 2 (-f: dj a 2Vlw - /^w?^) , 
^(pi)c b (P2)-c c (P3)c d ( P ,)) = -9 2 (-^jt cd -^fr bd -^r d dbc 

2 2 N 

_ a 4 robed _ ^± racbd _ a 4 radbc 

^(Pi)c\P2V(P.)c l {p,)) = - 9 2 {-^-ff d -^r M -^ff c 

<C>1)C^ 2 )C-^3)C^4)> = - 9 2 (j^tf* - ^fo^ - Yfo kM ) > 
(c a (pi)c fe (p 2 )c fe (p3) C '(P4)> = " <7 2 (" Yfo m - ^f^) , 

(c a (p 1 y(p 2 )c c ( P3 )c l (p 4 )) = o, 

(c\p 1 )c b ( P2 )c k (p 3 )c d ( Pi )) = 0, (A.2) 

where the momentum, pi, flow into each vertex and we have introduced the more compact 
notation, [T7) . 

jABCD _ pAB jiCD jABCD _ jeAB jeCD fj^ ^ 

for the four-point interactions with i, j, k and I denoting indices of objects in the group centre. 
The Feynman rules for the remaining possible combinations of diagonal and off-diagonal fields 
are trivially zero since they would generate group factors involving f^ c or which are both 
zero from the Lie algebra. 
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